Skyrmion Excitation in Two-Dimensional Spinor Bose-Einstein Condensate 
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We study the properties of coreless vortices(skyrmion) in spinor Bose-Einstein condensate. We 
find that this excitation is always energetically unstable, it always decays to an uniform spin texture. 
We obtain the skyrmion energy as a function of its size and position, a key quantity in understanding 
the decay process. We also point out that the decay rate of a skyrmion with high winding number 
will be slower. The interaction between skyrmions and other excitation modes are also discussed. 
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I. INTRODUCTION 

Topological objects have been attracting interest from 
various fields of physics for several decades 0. Roughly 
speaking, there are two kinds of topological excitations 
in two-dimensional space. The configuration of the first 
kind, such as vortex or monopole, has a natural singu- 
larity and depends on azimuthal angles even at infinity. 
These excitations have infinite kinetic energy unless they 
are coupled to other fields vanishing at infinity. In order 
to obtain a topological structure with finite energy, the 
configuration of the second kind must be uniform at in- 
finity, which means that the topological structure should 
be defined in a compactified space. The skyrmion is an 
example of the second kind topological structure, and the 
skyrmion excitation in n-dimensional space exists when 
the nth homotopy group of the internal space is nontriv- 
ial. 

Since its introduction in the 1960s in nuclear physics 2] 
and its application in QCD*^, skyrmions have been 
found in condensed matter systems such as quantum Hall 
effect and high temperature superconductivity. The 
achievement of Bose-Enistein condensate (BEc) in dilute 
Bose gases provide an opportunity to investigate many- 
body theory in new system. So far some topological 
excitations such as vortices and vortex rings in scalar 
BEc have been observed in a number of labs [3 ■ Recently 
the realization of a spinor BEc0, whose spin degrees of 
freedom are unfrozen has generated interest in studying 
richer topological excitations in such systems0|l]0]. 

The mean field description of spin-1 BEc was proposed 
in pioneer works of Ho and Ohmi et al 10]. The symmetry 
group of the ground state order parameter of an antifer- 
romagnetic spin-1 BEc is found to be C/ (1) x /Z^.pTj] 
where U{1) denotes the global phase angle and 5^ is a 
unit sphere denoting all orientation of the spin quanti- 
zation axis. The additional Z'^ is because the refection 
of spin quantization axis is equivalent to the change the 
global phase by tt. The general form of such a spinor field 
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where rft is the Bloch vector. It should be pointed out 
that the second homotopy group of 5*^ is homomorphic 
to integer group, this not only tells us of the existence of 
monopole excitation in three-dimensional antiferromag- 
netic BEc, which had been investigated by H.T.C. Stoof 
et al 8], but also implies the existence of the skyrmion 
excitation in two-dimensional case. 

It is naturally to ask if this excitation mode is ener- 
getically stable or not. In this paper, we answer the 
question within the Thomas-Fermi approximation and a 
variational method, and we find that the skyrmion is al- 
ways energetically unstable. In other words, in presence 
of any energy dissipation mechanism, the spin texture 
will always decay. However, this result does not imply 
that the skyrmion can not be created, a skyrmion with 
half topological charge has been successfully generated 
in a recent experiment by adiabatic deformation of the 
magnetic trap[i3|. Although the half skyrmion created 
in this experiment is different from that discussed in this 
paper, it is believed that a skyrmion with integer winding 
number can also be created in the near future. 

Furthermore, it is useful to discover whether the 
skyrmion decays by expanding to infinity size or shrink- 
ing to tiny size after it is created, and to understand how 
its center of mass moves. To answer these two questions, 
we need to obtain the energy of skyrmion as a function 
of its size and position, this energy function is the central 
result of our paper. We also find that the presence of a 
vortex influences skyrmion's motion. Additionally based 
on the results of Ref . 7i| , we will point out that the de- 
cay rate of high winding number skyrmions will be much 
slower than those with winding number 1. 



II. THE ENERGETIC STABILITY OF Q = 1 
SKYRMION 



When we only focus on the energy property of 
skyrmions, we can first neglect the interaction in spin 
channel which disrupt the S"^ order parameter, and the 
energy functional for such BEc can be simplified as fol- 
lowing: 
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H is the chemical potential and Vtrap is the confining trap 
potential. In this section, we will use conformal mapping 
to construct a general skyrmion excitation with winding 
number Q = 1. These variational wave functions are 
used to show that such skyrmions are always energetically 
unstable. 

Equation |5| is a nonlinear sigma model coupled to 
a (f^ model in an external potential. An n-skyrmion is 
an instanton solution to the 1 + 1 dimensional nonlin- 
ear sigma model|l|, this skyrmion minimizes the energy 
functional in the sector of each homotopy class, and can 
be constructed from the help of fractional linear mapping 
which maps the compactified complex space into itself, 
and the most general form of this mapping is 
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with the constraint aidi — biCi — 1, where N is the wind- 
ing number. Complex number uj are mapped to vectors 
on the unit sphere via 
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Owing to the conformal invariance of the nonlinear sigma 
model, its classical solutions are infinitely degenerate and 
the total energy is independent of the parameters in 
the analytical function / (z). The coupling between the 
spinor field ^ and the superfluid field ip breaks the confor- 
mal invariance, and skyrmions are not classical solutions 
to our model . The spin texture defined by equationiO 
can however be used in a variational calculation. The 
texture of spinor field contributes an effective potential, 
which changes the density profile. Given an effective po- 
tential, the minimum of the energy functional, is function 
of the parameters in the conformal mapping / (z). This 
function will tell us the information about the energetic 
stability of skyrmion excitation. 



We first consider the Q = I case, fl = f (z) ^ 
It is not difficult to show that 
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On the condition ad — be = 1, the effective potential is 
non-singular implying the density remains finite at the 
skyrmion core. This result is consistent with the coreless 
character of skyrmion excitation. The effective potential 
has two barrier localized at 2: = — - and z = —-, with 

c a ' 

the height |c| and \a\ respectively, where the density of 
the condensate should be smaller than the surrounding. 
Recall that in the repulsive case the interaction constant 
g is positive, both the interaction energy and the kinetic 
energy favor homogenous density profile, too much un- 
dulation will no doubt increase the energy, so the choice 
c — Q helps to decrease the energy. This value of c also 
produces the most symmetric texture. The function / (z) 
is reduced to a^(z — zi), where the parameters a and zi 
characterizes the size and location of the skyrmion. |l3j |. 

For simplicity, we first force zi = to be at the center 
of the trap. Non-zero zi will be discussed later. The effec- 
tive potential from the spin texture now is 7 — ^ — , 

[\i\' + \a\'rf' 

and the full potential is showed in Figure 1 for different a. 
Notice that the barrier height of the effective potential 
is proportional to |a|*, and that integrating the poten- 
tial energy over the whole two-dimension space results 
in a constant. In the limit a — > 0, the effective poten- 
tial spreads uniformly throughout the whole space. In 
the contrary limit of infinite a, the effective potential be- 
comes a S function. 

We now investigate the energetic stability of skyrmion 
in the framework of TF approximation where the term 
|V(y9p is neglected. In this approximation, the density 
profile is 
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in which ano = 

{^) and g = SttQscN, u! is the fre- The normalization condition and the constraint that the 
quency of the harmonic trap, and Q{x) is a step function. 
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condensate density is non-negative give the following two 
equations: 
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where x is the value of r at which the Thomas-Fermi den- 
sity vanishes. Solving the two equations we obtain the 
relationship between the chemical potential ^, Thomas- 
Fermi radius x and the skyrmion's size a. One can sub- 
stitute these relationships back to the energy functional 
and then obtain the minimal energy E as a function of 
size a which is plotted in Figure 2. 

Figure 2 shows that the minimal energy E has two 
minima occurring at zero and infinity, and there exists a 
critical size ac which corresponds to the maximal value of 
E. When a is quite large, the TF approximation fails, we 
use variational density profile to obtain a more accurate 
result [T3|. which is shown in Figure 3. 

Physically, we can understand the two figures in the 
following way. When the size parameter a is sufficiently 
small, the spin configuration within the TF radius be- 
comes uniform, the effective potential becomes fiat in- 
side the TF radius and just becomes a uniform shift of 
the chemical potential. So the energy approaches the 
ground state energy as a 0. In the opposite limit, the 
barrier is quite high and thus the density almost van- 
ishes at the center of the skyrmion, but at the same time 
the width of the barrier becomes narrow. The energy 
cost of the low density region is proportional to the vol- 
ume and decreases as a becoming larger. In this case the 
skyrmion can be observable by imaging the density pro- 
file, and observing the low density core. These figures tell 
us that the Q = 1 skyrmion is energetically unstable, in 
presence of any weak dissipation, it will either expand to 
become unobservable, or shrink to an infinitesimal size. 
As a semi-classical object the dissipative dynamics of a 
skyrmion depends on whether or not its initial size is 
larger than the critical size ac- 

For the same reason, the position zi of the skyrmion 
tends to move toward the edge of the condensate where 
the density is lower and the energy cost JVeff\ip\'^ is 
smaller. 



III. DISSIPATION DYNAMICS 

Although we believe that the energetically unstable 
skyrmion will decay in presence of any weak dissipation, 
the exact description of the dynamics depends on the de- 
tails of the dissipation mechanism. A possible source of 




FIG. 1: The trap potential added with the effective potential 
induced by Q = 1 skyrmion, in the unit of fuj, is plotted as 
a function of x/l for a/y'aHO = 1,1.5,2 respectively. Here 




FIG. 2: The minimal energy E in the unit of huj vs. the size 
of skyrmion a in the unit L = yfano- The parameter g/auo 
is set to 2. 



dissipation is interaction with the non-condensed com- 
ponent. Since skyrmion is a global topological object 
and the system is in a quite low temperature, this may 
not be important. Another source of dissipation is the 
spontaneous emission of other excitation modes, such as 
phonons mode or vortices, similar to the spontaneous 
emission of an excited atom. This mechanism is observed 
in numerical studies of superfluid vortex recconnection 



x/l 



FIG. 3: The minimal energy E in the unit of huj vs. size a in 
the unit of L, the result is obtained using variational method. 
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[l5j- Let us discuss the process in more detail. When a 
increases, the peak of effective potential grows and atoms 
are pushed away from the center of the skyrmion. The ki- 
netic energy of these atoms will lead to oscillatory motion 
of the condensate density. If the trap harmonic poten- 
tial is anisotropic, the extracted atoms can also rotate to 
form a vortex. These modes take energy away from the 
skyrmion, leading to a change of the skyrmion size. 



IV. INTERACTION BETWEEN A VORTEX 
AND A SKYRMION 

We rewrite (p as ^I^e*^. In equation Q the phase field 
9 is not directly coupled to the spinor field C, but they 
both interacts with the density field. Since both vortices 
and skyrmion have low density cores, it is favorable that 
the skyrmion core fits entirely within the vortex, thus one 
expects that the presence of vortex change the dominant 
trend of the evolvement of the skyrmion from expanding 
to shrinking, and one should find the skyrmion attracted 
to the vortex center instead of moving toward the cloud's 
edge. 

In principle, we can integrate out the ^E" field in the 
action and obtain the effective Hamiltonian to describe 
the interaction between a vortex and a skyrmion. Un- 
fortunately the 4'^ term makes this procedure difficult. 
However in the TF approximation, where the derivative 
terms of are neglected, the action is a quadratic form 
of the density field n. We can integrate it out and ob- 
tain an effective Hamiltonian describing the interaction 
between 9 and (: 
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The interaction energy will be smallest when the maxi- 
mum of \y9\ and |VC| coincide. This result is consistent 
with the above picture. 



with the height 
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this structure is markedly different from Vi whose peak is 
always localized at its center. In the center, the effective 
potentials is always zero for any n > 2. Figure 4 shows 
V2 for different value of a. As a is decreasing, the lo- 
cation of the off-center peak of the potential approaches 

/ 1 \ — ± 
infinity as (-) ^" and its height decreases as O" . Because 

the size of condensate is always finite, this implies that 
it is energetically favorable for a skyrmion to increase to 
a larger size so that the peak of the effective potential 
is outside the TF radius and the effective potential be- 
comes small and uniform inside the condensate. What is 
different from the Q = 1 case is the atoms must tunnel 
through the barrier of the effective potential as the size 
increases. Thus, the rate of the process will be charac- 
terized by the WKB tunnelling rate. The situation has 
been studied carefully in ReffT], and it showed that the 
tunnelling rate of such process may be long enough to 
form a dynamic mctastable skyrmion. 




FIG. 4: The trap potential added with the effective potential 
induced by Q = 2 skyrmion, in the unit of hui, is plotted as a 
function of x/l for a/^ano = 1, 3 and 5 respectively 



V. Q > 1 SKYRMION 



VI. CONCLUSION 



In this section we discuss skyrmions with Q > 1, 
demonstrating their differences from the Q = I skyrmion. 
Following the logic of the second section, we choose the 
most symmetric case, / (z) = a'^z", finding the effective 
potential 
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The barrier oi Vn [n > 2) lies on a circle 
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In summary, we constructed and studied the skyrmion 
excitation mode of the spinor field and showed that any 
finite size skyrmion is always unstable. Through a varia- 
tional study of the skyrmion excitation energy for differ- 
ent size and position, we find that it is energetically favor- 
able for skyrmion to expand to infinite size or shrinking to 
infinitesimal size, resolution in a uniform spin texture in 
most part of the condensate. We discussed the interplay 
between skyrmion modes and other excitation modes: (1) 
the emission phonons can dissipate the skyrmion energy 
and change its size, and (2) the vortex has attractive in- 
teractions with a skyrmion. We also showed that the be- 
havior of Q > 1 skyrmions are quite different from Q — 1 
skyrmion since the effective potentials induced by their 
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configuration have some markedly difference. However, 
all these discussions are concentrated on properties of 
the skyrmion mode its itself and a qualitative prediction 
on the skyrmion dynamics is made from the energetic 
consideration. To study the detail of skyrmion dynamic 
process, it is important to consider the interaction in spin 
channel, which cause spin fluctuations exceeding the S'^ 
internal space. ^] 

In the end, noticing that ( represents local spin-gauge 
degrees of freedom, we remark that the model we ex- 
plored in this paper shares some features of Yang-Mills 
theory, that is, the instanton solutions having the same 
winding number are saddle points of the Langrangian and 
are energetically degenerate although they have different 
shape and size, because of the local scaling invariance of 



the pure gauge theory. However, when the coupling to 
matter field breaks the conformal symmetry, the instan- 
ton solutions are no longer classical solutions and have 
different actions. In addition, the dynamic term of gauge 
field is absent in our model. The conclusion of energetic 
instability of any finite size skyrmion is determined by 
the above two features. 
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